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STABILITY OF TWO-FLUID WHEEL FLOWS 


by Eli Reshotko* and Carl F. Monnin 
Lewis Research Center 


SUMMARY 

The stability of an incompressible two-fluid wheel flow to infinitesimal 
helical disturbances is considered for the case where the inner fluid is heavy 
and the outer fluid is light. This situation may be viewed as a Rayleigh- 
Taylor problem in the frame of the rotating fluid and is dynamically unstable. 
As the density ratio of inner to outer fluid increases, the calculated growth 
rates increase and approach a limiting value as the density ratio becomes in- 
finite. Growth rates also increase with increasing axial as well as azimuthal 
wave number. Furthermore, the presence of a fixed boundary outside the light 
fluid tends to have negligible influence on the growth rates as the density 
ratio becomes very large. 


INTRODUCTION 

In some of the proposed gaseous nuclear rocket schemes, vortex containment 
is suggested for the maintenance of a critical mass of fissioning gaseous fuel 
with minimal losses (refs. 1 to 3). The present study is based upon Evvard’s 
concept of the wheel-flow reactor (ref. 3). In this reactor concept, a core of 
heavy fissioning gas in solid-body rotation is surrounded by an outer light gas 
coolant also in solid-body rotation at the same angular velocity. Further 
description of the hydrodynamic, nucleonic, and heat-transfer aspects of the 
wheel-flow reactor is given in reference 3. 

Unfortunately, the steady-state two-fluid wheel flow of a heavy gas core 
surrounded by a lighter, outer gas is unstable. The nature of this instability 
is the subject of this report. 

The stability of rotating flows is considered in detail by Chandrasekhar 
(ref. 4), but the reported treatments of wheel flow are limited to a single 
fluid in a bounded geometry. On the other hand, Chandrasekhar also considers 
the Rayleigh-Taylor stability problem for a stratified but nonrotating flow. 

Of interest here is the combination of these two cases which may in fact be 
viewed as a Rayleigh-Taylor problem relative to the rotating fluid. This fact 
was also pointed out by Melcher (ref. 5) in his general formulation for 
electrohydrodynamic waves in rotating systems. 

The wheel flow treated herein was idealized to facilitate obtaining solu- 
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Figure 1. - Effect of viscosity on growth rates of infinites- 
imal disturbances. 


Azimuthal 



tions; it consisted of two incompressible, immiscible fluids separated by a 
cylindrical interface. While the time -independent wheel flow is compatible 
with the complete Navier-Stokes equations including effects of viscosity, the 
viscous effects were omitted in this treatment of the stability of the config- 
uration. While this is a serious omission, the results yielded by this sim- 
plification are likely to be pessimistic. As pointed out by Chandrasekhar 
(ref. 4) for rotating fluids and also by Bellman and Pennington (ref. 6) in 
the case of the Rayleigh -Taylor problem, the effects of viscosity on the growth 
rates of infinitesimal disturbances tend to be as in figure 1. The salient 
effect of viscosity is a severe diminuation of the growth rate of short wave- 
length disturbances such that the maximum growth rate occurs at some finite 
axial wavelength. The results obtained herein should be reasonably correct in 
the long wavelength limit. The effect of azimuthal wave number on the growth 
rate will suggest whether large or small core fragments tend to break away more 
rapidly. 


BASIC FLOW 

The geometry of the two-fluid wheel flow is shown in figure 2. The problem 
is most conveniently handled in cylindrical coordinates. The axial or z- 
coordinate is directed into the paper. The inner fluid is given the subscript 
1, while the subscript 2 refers to the outer fluid. The interface between them 
is at radius a. The outer boundary at radius R bounds the outer fluid. 

The idealized basic flow is considered herein to be only azimuthal and 
solely a function of the radial coordinate. The following equations are the 
pertinent continuity and momentum equations, respectively: 
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Continuity: 



Radius, r* 

Figure 3. - Pressure distribution for inner 
fluid density greater than outer fluid density. 


sv* 

Te = 


o 


Radial momentum: 



_1_ 

p* dr* 


(la) 


(It) 


Azimuthal momentum: 


o = v *2£ + 1 


.* 


V 


\dr* 


?* dr* 


„*2 , 


(lc) 


(All symbols are defined in appendix A.) 


The general solution for azimuthal velocity 
obtained from equation (lc) is 

V*=Ar*+^ (2) 

r 


The wheel flow V = fir is the pertinent solution for the inner fluid. If the 
outer boundary is assumed to rotate at velocity J2R, then V = ftr is the com- 
plete solution and is compatible with the Navier-Stokes equations including 
viscosity. The pressure distribution in the wheel flow from solution of equa- 
tion (lb) is 


P = 


* , 

P n + Pi 


2 *2 
Q r 


P* + P s!s! + p aHr * 2 - a 2 I 
*0 H 1 2 2 2 


0 < r* < a 
a < r* < R 


(3a) 

(3b) 


For a heavy fluid core surrounded by a light fluid, the inertial nature of the 
instability is apparent from the pressure distribution for > p 2 (fig. 3). 

If an element of heavy fluid from just inside the interface (r* = r£) is dis- 
placed to just outside the interface (r* = rf) while conserving its angular 

momentum, a pressure gradient dp*/dr*= p^/r 2 [y^(r*/r 2 )J ( represented by the 
dashed line in fig. 3) is required for radial balance. However, the available 

pressure gradient at r| is dp*/dr* = P 2 (V g /r 2 ) (solid line), which for 
wheel flow is p 2 /r 2 [v^(r 2 /r*)J . Thus in the limit of infinitesimal displace- 
ment [(r 2 - r* ) /r"£ « lj as long as p^/p 2 > 1 the heavy fluid encounters a 
pressure gradient inadequate for radial balance on penetrating the interface 
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and the element of heavy fluid will continue its motion outward. This behavior 
is characteristic of Rayleigh -Taylor instabilities 


FORMULATION OF STABILITY PROBLEM 


Disturbance Equations 

The equations governing infinitesimal disturbances are obtained from a 
linearization of the complete equations of motion about the assumed basic flow. 
As previously mentioned, only inviscid disturbance motions will be considered 
in this treatment. For the basic wheel flow V* = ftr*, the disturbance equa- 
tions for each of the fluids are: 


Continuity: 


Momentum : 


1 5 / # -x- \ , 1 5v~* , 5w~* n 

— — ( r *u* ) + — + — = 0 


Sr 




^ + a 2s* . 2Uv* 


St 


* 


be 


Sv* dv* * 

+ zr + 2flu 


1 bx*_ 
~ P Sr* 

1 S?t* 

pr* 


Sw* Sw* _ 1 Sfl* 

at* de ' ' » a 2 * 


(4a) 


(4b) 

(4c) 

(44) 


These equations may also be written in dimensionless form for later convenience. 
When the disturbance quantities and the independent variables are referred to 
the following reference values , 


v* ef = fla 

(5a) 

I^ref = a 

(5b) 

* _ ^ 0 2 q 2 
P re f a 

(5c) 


^Ehe argument presented here parallels that given by von Karman (see Lin, 
ref. 7) in support of Rayleigh's stability criterion for rotating fluids. The 
equivalent to Rayleigh's criterion for a stratified incompressible ideal fluid 
(inviscid and not thermally conducting) is that the flow configuration is 
stable if d/dr(pV 2 r2) > o. The two-fluid wheel flow under consideration 
(p n /p_ > l) is clearly unstable by this criterion. 

Ci 
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t* - i 

^ref - o 


(5d) 

^ef = ^ 


(5e) 

The following equations are the dimensionless disturbance 

equations : 


1 S / \ , 1 dv dw n 

r 5? (ru) + T Se + 5J ■ 0 


(6a) 

3u , du ~ 

5t + Se - 2v ' - 5? 


(6b) 

bv , bv , „ 1 Sit 

5t + Se + 2u = -?S? 


(6c) 

soVo 

II 

+ 


( 6d) 

Since the coefficients of equations (6) are independent of t, 0, and z, 

the equations allow disturbances of the form q(r )e i ^ kz+m0 " C0t ^ , where q(r) is 
a complex disturbance amplitude, k and m are the axial and azimuthal wave 
numbers, respectively, and is a complex frequency. The wave numbers k 

and m are real numbers; m is an integer. The real part of cd is the rota- 
tional frequency of the disturbance flow, while the imaginary part <£>± repre- 
sents the amplification rate. Disturbances are growing, neutral, or decaying 
according to whether is positive, zero, or negative, respectively. For 

the assumed disturbance form, equations (6) become 

u' + u + i2E +ikv= o 
r r 


(7a) 

icuu + 2v = it T 


(7b) 

^ imrt 

icov - 2u = — 


(7c) 

io)w = ikrt 


( 7d) 


where primes denote differentiation with respect to r and £5 = a> - m is 
related to the angular velocity of the disturbance pattern relative to the mean 
rotation.^ Thus, to an observer at a given axial location z and moving with 
the mean rotational velocity ft, the problem may be viewed as a Rayleigh-Taylor 

2 In terms of dimensional quantities, <x£/m = o/ t* - ft where (0/t*) is 

the absolute angular velocity of the disturbance pattern. To an observer at a 
given axial station z, the quantity = (m/ft)( $/t* - ft) is proportional to the 

angular velocity of the disturbance pattern relative to the mean rotation. The 
dimensionless angular velocity of the wave is simply e/ftt = 1 + oS^/m. 
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problem with a centrifugal rather than a gravitational driving force. 

From equations (7b), (7c), and (7d), the velocity fluctuation amplitudes 
may be expressed in terms of jt : 


(8a) 

4 - 



a*- - as 


V = 


4 - co 


(8b) 


w = 




CO 


(8c) 


Substitution of equations (8) into the disturbance continuity equation (7a) 
yields the differential equation for the pressure fluctuation amplitude: 


♦ r - H + - 


0 


(9) 


where 


A 2 = 1 - — 

err 


( 10 ) 


The general solution to equation (9) is any linear combination of Bessel 
functions of the first and second kind of the proper order and argument. The 
form used herein is 


n = AJ m (ikAr) + BH^ 1 ) (ikAr) (ll) 

where the first term is regular at the origin while the second is regular at 
infinity, provided that the real part of kA is positive for large values of 
kA. 


In the inner fluid, the second term is inadmissable so that 

Jtj_ = A 1 J m (ikAr) (12) 

In the outer fluid, the general solution is 

*2 = A 2 J m (ikAr) + B 2 H i l} (ikAr) (13) 

If the outer boundary is at infinity, the coefficient Ag must vanish as a 
consequence of the proper satisfaction of boundary conditions. 
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Boundary Conditions 


The boundary condition at r = 0 has already been applied. The remaining 
conditions are that the radial velocity must vanish at r = R and that the 
radial velocity and normal stress be continuous at the interface r = 1 
(r* = a). These conditions are treated consecutively. 

Radial velocity u = 0 at r = R . - Pbrom equations (8a) and (13), 



where 


and 



SmH^ 1 ^ ik£r) 

(Ur 


(15) 


(16) 


As R approaches infinity., approaches zero, as seen from equa- 

tion (14). 

Radial velocity u^ = u 2 at r = 1. - From equations (8a), (12), (13), 
(15), and (16"), this condition yields the following relation between A ^ and 
b 2 : 

A 1 \(B) 

B / (l) ■ 7 (R) U j 

2 "m * m 

Continui ty of normal force at interface . - The satisfaction of this 
boundary condition is accomplished by obtaining the sum of pressure plus cen- 
trifugal forces for each fluid at the displaced location of the interface and 
then equating these sums. Because two separate fluids with different densities 
are involved, the force balance is carried out in physical variables: 
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Po 


G 2 a 2 


t 1 + s l (1) ] 


+ p^fl 2 a 2 3t^(l) 


= P 0 + Pi + ~ 2 5 jt 1 + & 2 (1) ] 2 “ X | + P 2 « 2 a 2 rt 2 ( 1) (18) 

When squares of disturbance quantities are dropped, the linearized form of 
equation (18) becomes 

p l[^l (l) + *1 (1 )] = + * 2 d)] (19) 

The amplitude of the radial interface displacement % is obtained by a time 
integral of the velocity in the rotating frame; that is. 



( 20 ) 


General Dispersion Relation 

After the appropriate substitutions are made into equation (19), the 
following dispersion relation is obtained: 


V ikA) & m (R) J m ( ikA) 

</m^ 1 "l R ^ 1 ^ P 1 J m^ 

SjR) ‘ P 2 A (lj 

J 7 ^ t 


( 21 ) 


where the 34 m and 7m functions are those defined in equations (15) and (16) 
and A is defined in equation (10). The dispersion relation (21) is a complex 
equation. In fact, since the complex unknown 26 appears in the .arguments of 
the Bessel functions through A, solutions to the dispersion relation (21) are 
probably best obtained by an iterative trial-and-error procedure. Such a 
procedure has been employed to obtain numerical solutions for arbitrary complex 
frequencies where the Bessel functions for arbitrary complex arguments are 
evaluated by the method of reference 8 . 


The dispersion relation (21) has multiple roots, not all of which have 
physical significance. This difficulty arises particularly because of the 
assumption regarding the positiveness of the real part of A that was made in 

choosing the Hankel function of the first kind rather than the Hankel 

function of the second kind H^ 2 ) for the general solution (ll). The appro- 
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priate roots of the dispersion relation will he sorted out hy examination of 
various limiting solutions of the dispersion relation (2l). 

Before proceeding, some remarks are in order on the effect of the outer- 
most boundary at r = R. The boundary condition that u = 0 at r = R enters 
the dispersion relation (2l) through the terms in the first group- 

ing of the denominator of the right side. For large inner to outer-density 
ratio p 2 _/Pg » the first portion of the denominator is small compared with 

that term which has p-j/pg as its coefficient, and its effect on the result 
is very weak. Furthermore for large arguments (k£R), CHn(R)// m (R)] ~ e- 2k ^ 
becomes quite small. For the values of p^/pg and R of interest in wheel- 

flow reactors, the stability charactersitics are essentially those of the un- 
bounded configuration. The solutions presented herein are, therefore, for the 
unbounded configuration R -> oo. 


SOLUTIONS TO DISPERSION RELATION FOR UNBOUNDED CONFIGURATION 


For the unbounded configuration, 
relation (2l) reduces to 


[3*m(R)i/m(R)] 


0 and the dispersion 


2£ - 4 = 




p l 


1 

p 2 


fl^ 1 ) (ikAr)l 

L m 'J r=1 2m 

J m (UcAr) 

1 

r=l 2m 


(ikA) 2 J m ( ikA) 2 


( 22 ) 


Solutions to dispersion relation (22) are obtained and discussed for 
axially symmetric disturbances (m = 0) and then for the helical disturbances 
with flute number m different from zero. 


Axially Symmetric Disturbances 

For axially symmetric disturbances (m = 0), equation (22) becomes 




iH^ 1 ^ (IkA) 

+ h 

J 0 (ikA) 

-H^ 1 ^ (ikA) 

P 2 

-U-lC IkA) 


(23) 
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In this particular case 2S = oo since m = 0. Solutions are now sought subject 
to the condition that the real part of A is greater than zero. Assume for 
the moment that for m = 0, A is real and positive. Then, the bracketed 
ratios of Bessel functions in equation (23) are real and positive so that the 
right side of equation (23) is real but always negative. Since the left side 
of equation (23) is just o?A 2 , a solution is possible for negative real values 
of or purely imaginary values of go, which are obtained for 


^(1) (ikA) 


kA > 


-Hp^MikA) 


J 0 (ikA) 
L - !!-,_( ikA)_ 


Pi 

— - 1 

p 2 


(24) 


From the definition of A 2 , these are the only admissible solutions with 
A r > 0. Thus the axially symmetric disturbances may be described as standing 
waves (o^ = 0) that are amplified (coj_ > 0).^ 


• Solutions to dispersion relation (23) were obtained (by hand calculation) 
for density ratios p^/pg of 2, 10, 100, and oo. These results are shown in 
figure 4. As expected from the qualitative discussion in the INTRODUCTION, the 
growth rates increase with increasing density ratio. For density ratios p^/pg 
greater than about 10, the results are not very different from those for an 
infinite density ratio, which indicates that at these density ratios the pres- 
sure in the outer fluid has negligible effect in decelerating a displaced fluid 
element. For large axial wave numbers , the results are well represented by the 
asymptotic solution to equation (23) for large kA obtained in appendix B, 
namely. 



(25) 


(26) 


(27) 


3 Both amplification (u>j_ > 0) and decay (ay. < 0) of the disturbances are 
possible. Only the possible amplification, however, is relevant to the assess- 
ment of the stability of the configuration. 
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Solutions for very long wave- 
lengths k -» 0 are obtained at 
finite kA at the value corre- 
sponding to the equal sign in equa- 
tion (24). Since kA for axially 
symmetric disturbances is always 
greater than 4, the solutions (25) 
to (27) are good representations of 
the results of figure 4 for all 
wave numbers and density ratios • 


Nonaxially Symmetric 
Disturbances 


For the nonaxially symmetric 
disturbances, the dispersion rela- 
tion (22) is solved numerically by 
trial and error on an IBM 7094 com- 
puter with the pertinent Bessel functions for arbitrary complex arguments cal- 
culated according to the method of reference 8. The results for density ratios 
2, 10, and 100 are shown in figure 5. The numerical results of these densities 
agree well with asymptotic solutions to dispersion relation (22) in both the 
long- wave length k -> 0 and the short -wave length limits (appendixes C and D, 
respectively) . 

Limiting solution for long wavelength * - Long-wavelength solutions k -> 0 
for nonaxially symmetric disturbances are sought in the limit kA -» 0, despite 
the aforementioned result for m = 0. The pertinent solution (from appendix C) 
is 


Ratio of inner to 
outer fluid 
density, 

Nh 



0 2 4 6 8 10 


Axial wave number, k 

Figure 4. - Growth rates of axially symmetric disturbances in un- 
bounded configuration. 


“r “ ' £ 


(28a) 


or 


^ = m - & 


(28b) 


and 


“j. = 


■y/mSt - 1 

£ 


(29) 


where £ is defined in equation (27). This solution satisfies the condition 
A r > 0 ( appendix C ) . 

Limiting solution for short wavelength . - From appendix D, the leading 
terms of the solution of the dispersion relation (22) in the limit kA -» oo are 






where 


( 4 + i^ + • • •) (s3) 

In these expressions, kA is a real number. Note that o >^ q is independent of 

the flute number m and is identical to the asymptotic result for axially 
symmetric disturbances. The dependence of (&± on flute number m is of higher 
order , as is the departure of kA from the real axis. 

The numerical solutions presented in figure 5 are entirely consistent with 
these limiting forms. The nonaxially symmetric disturbances are growing waves 
that for long wavelengths rotate at 1 - l/m£ of the angular velocity of the 
basic flow. As the wave number increases , the angular velocity of the wave in- 
creases asymptotically and approaches the basic flow angular velocity in the 
short-wavelength limit. For kink instabilities (m = l), the increase in wave 
angular velocity occurs rather abruptly at dimensionless axial wave numbers 
between 0 and 1. The waves at higher flute numbers rotate at higher fractions 
of basic flow angular velocity and approach that velocity more gradually with an 
increase in wave number. The growth rates themselves generally increase with an 
increase in flute number m at all axial wave numbers. An exception occurs 
for m = 1 near k = l/2 for the density ratios 10 and 100. This exception 
is as yet unexplained. As expected, the level of growth rates increases with 
increasing density ratio in about the same way that it did for axially 
symmetric disturbances. For large axial wave numbers, the growth rates ap- 
proach the cd^o given by equation (33) regardless of flute number. 


RELEVANCE AND SIGNIFICANCE OF RESULTS TO WREEL-FLOW REACTOR 

The wheel-flow reactor (ref. 3) consists of rotating hot compressible 
gases, where the inner fissionable fluid may be quite a bit hotter than the 
outer light gas propellant. The interface between these two fluids will not be 
sharp because of the diffusion of heat and because of species concentration 
diffusion. The treatment of such a flow in terms of two constant-density 
immiscible fluids is sufficiently idealized that only qualitative statements 
may be made to relate the results of the present stability calculations to an 
actual wheel-flow reactor. Such qualitative statements are expected to have 
some validity because, as previously stated, the instability is primarily iner- 
tial. The statements in this section are based on the premise that the normal 
mode disturbances with the largest growth rates are expected to be the most 
troublesome. 


For conceivable molecular weights and 
ratios are in the range 10 < P^/pg ^ 1°G. 

tive difference exists between the results 
of density ratios. Significant reductions 

p /p <2. 

l' 2 ~ 


temperatures, the pertinent density 
Little qualitative or even quantita- 

(figs. 5(a) and (b)) in this range 
in growth rate are obtained only for 
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Although calculated growth rates increase with increasing axial wave 
number in the present inviscid calculations, the large wave-number disturbances 
(short wavelength) will be damped by the action of viscosity. This has been 
verified by the present authors in a preliminary calculation not presented 
herein. 

Accordingly, the disturbances most likely to be amplified are those with 
axial wavelengths larger than the interface radius, a (k < 2 jt) . For this range 
of axial wave numbers, the growth times are of the order of the period of basic 
flow rotation and increase with azimuthal wave number or flute number. This 
increase indicates a rather rapid growth of disturbances , with smaller core 
fragments tending to break away most rapidly. Same of these troublesome long- 
wavelength disturbances can perhaps by eliminated by choosing the reactor 
length to be short, of order jta, so that the minimum axial wave number is 
k = 1. 

The previous statements should be reexamined with the aid of a less ideal- 
ized study that might include effects of compressibility and species diffusion. 
The fissionable core would be substantially ionized and electrically conducting. 
This suggests the possibility that additional stabilization might be obtained 
by use of imposed axial magnetic fields. 

SUMMARY OF RESULTS 

The stability of an incompressible two-fluid wheel flow to infinitesimal 
inviscid helical disturbances has been considered for the case where the inner 
fluid is heavy and the outer fluid is light. The following are the important 
results : 


1. For large inner- to outer-fluid density ratio, the presence of a physi- 
cal boundary outside the light fluid has only a very weak effect on the stabil- 
ity characteristics. For the values of density ratio and outer radius of in- 
terest in wheel-flow reactors, the stability characteristics are essentially 
those of a configuration where the outer fluid is unbounded; this applies both 
to the axially symmetric disturbances and to the nonaxially symmetric disturb- 
ances. 


2. The axially symmetric disturbances are identified as standing waves 
that are amplified. The growth rates increase with increasing density ratio. 
For density ratios greater than about 10, the results are not very different 
from those for infinite density ratio. 

3. The nonaxially symmetric disturbances are helically propagating waves 
that are also amplified. For a given density ratio, the growth rates increase 
with azimuthal wave number. The magnitude of the growth rates varies with den- 
sity ratio in about the same way as for axially symmetric disturbances. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, May 24, 1965 . 
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APPENDIX A 


A,B 

a 

34 

/ 

k 

L 

m 

CM 

p 

q.(r) 

R 

£ 


r 

t 

u 

V 
v 
w 
z 

A 

e 

V 


SYMBOLS 


constants in solution of differential equations 

radius of two-fluid interface 

defined by eq. (15) 

defined by eq.. (16) 

axial wave number 

reference length 

azimuthal wave number or flute number 

order of magnitude of x 

pressure 

complex disturbance amplitude 
radius of outer boundary 



radial coordinate 
time 

radial disturbance velocity 
azimuthal velocity of basic flow, Or 
azimuthal disturbance velocity 
axial disturbance velocity 
axial coordinate 

defined by eq. ( 10 ), A 2 = 1 - -77 

OCr' 

azimuthal coordinate 
kinematic viscosity 


| radial displacement of disturbance flow 

jt disturbance pressure 

p density 

SI angular velocity of solid body rotation 

co complex angular disturbance frequency 

co co - m 

Subscripts : 

i imaginary part 

r real part 

ref reference quantity 

1 inner fluid 

2 outer fluid 
Superscripts : 

dimensional quantity 

1 differentiation with respect to r 
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APPENDIX B 


SOLUTION FOR LARGE WAVE NUMBER OF DISPERSION RELATION FOR AXIALLY 

SYMMETRIC DISTURBANCES IN UNBOUNDED CONFIGURATION 

The dispersion relation pertinent to axially symmetric disturbances in the 
unbounded configuration is that of equation (23), namely, 


or 


4 


- 1 (kA) 


iH^ ) 1) ( ikA)~ 

+ ei 

J 0 (ikA) 

-Hj 1 ) (ikA) 

P2 

-iJ x (ikA) 


(23) 


For large values of kA, the functions in brackets are approximately 


iH^ l} (ikA) , 

Y . - 1 _ + 

r (l) 2kA 


-Hi 


(ikA) 


(Bl) 


j q( ikA) ! 4 _ 

- iJ 1 ( ikA) 2kA 


(B2) 


The dispersion relation (23) may be written 


Qj — 4 — ~ 


kA 

IT 


( 1 " 2lkA + ‘ * •) 


(B3) 


-where 


ft - 



(27) 


Since, in this case, kA is real and is purely imaginary, 



Now from the definition of A (eq. 10), 


(25) 
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Thus, 


The wave 


~ 2.2 
CD A = 


£j A 

CD - 4 


(B4) 



number k is then given by the expression 



(B5) 


(26) 
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APPENDIX G 


SOLUTION OP DISPERSION RELATION FOR UNBOUNDED 
CONFIGURATION FOR SMALL WAVE NUMBER 

The dispersion relation for the unbounded configuration R -* oo (eq.. (22)) 
is to be solved in the limit of small wave number k. The first step in so 
doing is the evaluation of 

J m ( ikA) 

and 


(ikAr) 
H^ 1} (ikA) 


r=l 


for small values of the argument kA. When the appropriate relations are used 
for kA -* 0, 


[ J m (^)]) =1 

J m (lkA) 


= m + 


(Cl) 


and 

(ikAr) 

L m Jr=l 

■ • — / - \ — = -m + . . . 

< 1) (ikA) 

The dispersion relation (22) for small kA may be written 

(a? - + 2a> + m) = 0 


where 


(C2) 


(C3) 


R s 



(27) 
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The solutions of this dispersion relation are 


CD = ± 2 


(C4) 



(C5) 


For solutions (C4), A 2 = 0. These solutions are rejected since the condi- 
tions A r > 0 is not satisfied, namely, that disturbance amplitudes must decay 
exponentially as the radius becomes infinite. 


Solutions (C5 ) are both acceptable; that with a positive imaginary part 
indicates growth of the disturbance, while the other indicates decay. Both are 
possible, but only the growth solution is of interest here since it persists. 
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APPENDIX D 


SOLUTION OF DISPERSION RELATION FOR UNBOUNDED 
CONFIGURATION IN LARGE WAVE NUMBER LIMIT 


For large values of the argument kA, 

[j m (ikAr)]^ i 




and 


fc l) ( ikAr )1 
L -*r=l 

H (ikA) 
m 


( 1 -2S + ' ' •) 


= - kA ( 1 + is + • • •) 


(Dl) 


(D2) 


Equations (Dl) and (D2) are independent of the flute number m, to the order 
considered, so that the dispersion relation (22) may be approximated as 


00 



(D3) 


The last term in (D3) is of order l/(kA)^ compared with the leading term on 
the right side of (D3)> and it will be neglected. 

Equation (D3) may now be alternatively written 


~ A , ~ V Cl ill 

“ “ “ 1,0 = ^2 


2 

where q represents the grouping of terms 

2 kA / 

“i,o - W ~ ( 4 + 


(D4) 


(D5) 


The right side of equation (D4) is always of unit order. For large values 
of kA (therefore, large values of q), the left side becomes very large and 
suggests a solution of the form 

So = ioo-j^Q + ^1 (D6a) 
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where 


a^L 

%0 


« i 


(D6b) 


Substituting (D6a) into (D4) and keeping only the leading terms of the expan- 
sion yields 


= - 


m 




(D7) 


which is seen to satisfy (D6b) for kA » 1 and is the leading contribution to 
the real part of 25 . 

To obtain the asymptotic form of the wave number k, the quantity A must 
first be evaluated. From its definition. 



( 10 ) 


With the results (D6a) and (D7), 



(D8) 


When terms are kept up to order 



compared with the leading terms 


0 ) = 10 ). 


m 


1,0 ’ * 2 4,o 


(D9) 




A = /l + 


CJD? 


i,0 


(DIO) 


and 


k = 


kA 


r 


4,o 


(Dll) 


To this order of approximation, A is a real quantity, and, therefore, kA 
is real as well. 


24 



REFERENCES 


1. Kerrebrock, Jack L.; and Meghreblian, Robert V.: Vortex Containment for the 

Gaseous -Fission Rocket. J. Aerospace Sci., vol. 28, no. 9, Sept. 1961, 
pp. 710-724. 

2. Rosenzweig, M. L.; Lewellen, W. S.; and Kerrebrock, J. L. : The Feasibility 

of Turbulent Vortex Containment in the Gaseous Fission Rocket. ARS J., 
vol. 31, no. 7, July 1961, pp. 1873-1883. 

3. Eward, John C.: Wheel-Flow Gaseous-Core Reactor Concept. RASA TN L-2951, 

1965. 

4. Chandrasekhar, S.: Hydrodynamic and Hydromagnetic Stability. Oxford 

Univ. Press, 1961. 

5. Melcher, J. R. : Electrohydrodynamic Waves in Rotational Systems. Quart. 

Prog. Rept. No. 69, Res. Lab. of Electronics, M.I.T., Apr. 15, 1963, 
pp. 102-113. 

6. Bellman, R.; and Pennington, R. H. : Effects of Surface Tension and 

Viscosity on Taylor Instability. Quart. Appl. Math., vol. 12, 1954, 
pp. 151-162. 

7. Lin, C. C. : The Theory of Hydrodynamic Stability. Cambridge Univ. Press, 

1955, pp. 49-50. 

8. White, J. E. , Jr.: A FORTRAN Program to Calculate Bessel Functions of 

Integral Index and Complex Argument. Rept. No. R-59, Brown Eng. Co., Inc., 
July 1963. 


NASA-Langley, 1965 E — 2936 


25 



"The aero?2autical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space . The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 
important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless 

of importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distri- 

bution because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Technical information generated in con- 

nection with a NASA contract or grant and released under NASA auspices. 

TECHNICAL TRANSLATIONS: Information published in a foreign 

language considered to merit NASA distribution in English. 

TECHNICAL REPRINTS: Information derived from NASA activities 

and initially published in the form of journal articles. 

SPECIAL PUBLICATIONS: Information derived from or of value to 

NASA activities but not necessarily reporting the results of individual 
NASA-programmed scientific efforts. Publications include conference 
proceedings, monographs, data compilations, handbooks, sourcebooks, 
and special bibliographies. 


Details on the availability of these publications may be obtained from: 


SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



